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ABSTRACT 


Dynamic  scaling  theory  is  used  to  obtain  expressions  for  the  far  infrared  absorption  coefficient 


a  of  metal-insulator  granular  composite  materials  near  the  dc  percolation  threshold  pc  in  the 

A  •  • :  '  -  •• 

quasistatic  limit.  Anomalous  frequency  dependence  is  predicted  near  p^  for  50  iT  radius  Au  particles. 
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The  frequency  dependence  of  -a'  obtained  from  an  analysis  of  published  data  on  Au  blacks  is 
inconsistent  with  the  scaling  theory  for  Drude  Au  particles  in  the  three  dimensional  limit,  and  the 
theoretical  estimate  of  the  magnitude  of  a  at  w  =  100  cm-1  is  a  factor  of  seven  too  large. 
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The  anomalous  enhancement  of  the  measured  far  infrared  absorption  coefficient  of  granular 
metal-insulator  composite  materials  of  low  metallic  content  with  respect  to  theoretical  predictions  has 
attracted  interest  for  over  a  decade1-'.  After  clustering  was  convincingly  associated  with  enhanced 
absorption  experimentally4-6,  theoretical  models8-10  have  focused  on  mechanisms  based  on  clusters  of 
particles.  Some  of  the  proposed  mechanisms  1  treat  tenuous  clusters  and/or  clusters  with  metallic 
volume  fractions  very  close  to  the  dc  percolation  threshold  pc. 

Effective  medium  theories  for  metal-insulator  composite  materials  fail  in  the  region  of 
composition  near  pc.  The  scaling  theory  of  percolation11,  based  on  the  universality  of  behavior  near 
pc,  has  been  generalized  to  treat  the  frequency-dependent  complex  effective  dielectric  function12-14. 
Recently,  Niklasson  and  Granqvist15  (NG)  argued  that  scaling  theory  applies  to  the  mid-  and  far- 
infrared  absorption  coefficient  extracted  from  data  on  evaporated  Au  blacks16-18. 

In  this  paper,  we  apply  dynamic  scaling  theory  to  obtain  expressions  for  the  far  infrared 
absorption  coefficient  a  of  a  granular  metal-insulator  composite  material  for  p  near  pc.  We  predict 
anomalous  frequency  dependence  very  close  to  pe.  We  examine  the  possibility  of  observing  the 

O 

predicted  effects  experimentally  in  the  far  infrared  for  50  A  radius  Drudc  Au  particles  in  a  perfect 
insulator.  It  is  important  to  verify  that  the  quasistatic  approximation  is  valid.  Finally,  we  reexamine 
Niklasson  and  Granqvist’s15  analysis  of  the  data  of  Harris  et  a/.16-18  and  find  that  our  expression  is 
inconsistent  in  frequency  dependence  if  the  sample  is  in  the  three-dimensional  (3D)  limit.  Also,  the 
predicted  magnitude  of  a  exceeds  the  measured  value  by  a  factor  of  seven.  We  conclude  that  further 
experimental  work  on  well-characterized  samples,  such  as  cermet  films,  is  required  to  test  the  scaling 
theory. 

Consider  small  metallic  particles  with  complex  dielectric  function  Crn(w)  imbedded  in  a 

nonabsorbing  host  with  dielectric  constant  Cq.  If  (p  —  pcj  =  |Ap|<C  1  and  Ifp/tml*^^  the  singular  part 

13  19 

of  the  effective  complex  dielectric  function  c(w;p)  can  be  written  in  the  scaling  form 
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i/cm  =  |Ap|l  F±gjJ  lApp8-1),  (1) 

where  the  subscripts  +  and  —  refer  to  p>pc  and  p<pC)  respectively.  Define  z  =  (fp/cm)  |Ap|  — s  —  l. 
Following  Bergman  and  Stroud13,  for  |z|4Cl,  F_^(z)  can  be  expanded  to  obtain 


e  —  B—  Cq  |Ap|  s  +  C—  (f^Am)  |Ap| 


—  t  — 2s 


for  p<pc  and 

c  =  A+  cm  | Ap +  B_j_  fj-j  |Ap|  s  +  ...  (3) 

for  p>pc.  The  coefficients  A^_,  B and  C—  are  presumed  real  since  the  theory  should  apply  when  em 
and  Cq  are  real.  Very  close  to  pc,  |z|»l.  Since  ?(w;p)  is  independent  of  p  to  lowest  order, 

t/(s+t)  s/(s+t)  ... 


The  absorption  coefficient  o(w;p)  is  defined  by 


o  =  2  (w/c)  Im{4?}, 


where  c  is  the  speed  of  light.  Define  £  =  ?j^+i?j.  Let  fm(w)  be  described  by  the  Drude  model,  which 
applies  to  many  metals  in  the  far  infrared.  If  wr<C  1,  cm  ~  iwp r/w,  where  r  is  the  relaxation  time 
and  Wp  the  plasma  frequency.  For  |z|<£l  with  p<pc,  since  £j^~z  and  fj~z2.  Then 

o~(w/c)  «j/  so 
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To  lowest  order,  a  is  independent  of  p,  but  shows  anomalous  frequency  dependence.  Equations  6  —  8 
are  valid  if  the  singular  part  of  c  is  the  dominant  contribution. 

To  examine  whether  the  predicted  effects  are  observable  in  the  far  infrared,  we  consider  a=50 
A  radius  Drude  Au  particles  supported  in  vacuum  (fp  =  l).  For  Drude  Au,  the  plasma  frequency20  is 
<Lp  =  7.27xl04  cm-1  (w  =  w/27tc  is  the  frequency  in  wavenumbers.).  For  very  small  particles,  the 
carrier  scattering  is  dominated  by  the  surface  of  the  particle21,  so  we  make  the  approximation 
r~a/vp,  where  Vp  is  the  Fermi  velocity.  Thus,  l/(2ircr)  =  1490  cm-1  (Ar  =  6.73 /rm).  wr<l  is 
reasonably  obeyed  for  w<100  cm-1.  At  w  =  100  cm-1,  |ej-j/crn|  =  2.82xl0~5 <g;  1 ,  as  required. 

The  quasistatic  approximation,  which  underlies  this  scaling  theory,  applies  if  k£<Cl,  where 
k  =  2x/A  =  27ru )  (k  is  the  wavevector,  A  the  wavelength)  and  £  is  the  correlation  length  from  percolation 
theory22.  Clearly  £^A  sufficiently  close  to  pc,  so  it  is  important  to  verify  that  one  can  achieve  the 
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conditions  required  to  observe  anomalous  absorption  within  the  bounds  imposed  by  the  quasistatic 
limit.  From  percolation  theory22,  £~a  j Ap|  v .  We  have  introduced  the  particle  radius  a  as  a  length 
scale  for  numerical  estimates.  In  the  quasistatic  limit,  (2eC(2tr£)  1  =|Ap|t,/(27ra)  =  Wy  The  crossover 
frequency  w*  for  anomalous  frequency  dependence  is  obtained  by  setting  |z|=l.  The  result  is 


/  =  ^  |AP|‘+*. 


is  required  for  anomalous  frequency  dependence.  Eqs.  6  and  7  apply  when  More 

restrictive  bounds  on  the  region  of  anomalous  frequency  dependence  are  obtained  by  considering  the 
wavelength  inside  a  cluster  of  size  £,  Ac  =  A/nc,  where  nc  =  Re|^tc|  and  cc  is  the  dielectric  function  of 
the  cluster.  Although,  the  quasistatic  approximation  is  easily  satisfied  for  an  isolated  Au  particle,  if  we 
assume  that  the  cluster  is  a  Au  sphere,  we  obtain  the  very  strong  requirement  u<g.u>0,  where 


“2  =  f  (s%)2  lAP|2l/- 


It  is  more  plausible  to  approximate  the  dielectric  function  of  the  cluster  by  the  effective  dielectric 
function  at  pc  (Eq.  4).  The  bound  becomes  where 


w3  =  [(a/c)cos(|  ^)] 


2(s+t)  —  t  — s  2i/(s+t) 

s+2t  s+2t  ,  2  \sH-2t  |  a  |  s+2t 
CD  (wpr)  |Ap| 


Figure  1  shows  the  various  regions  of  behavior  near  pc  for  50  A  radius  Drude  Au  particles  for 
both  3D  and  2D.  Following  NG,  for  3D  we  use23  s=0.75  and24  t  =  1.94.  Also25,  r/  =  0.88.  For  2D  we 
use23  s  =  t  =  1.3  and26  u  =  4/3.  We  conclude  that  the  possibility  of  observing  anomalous  frequency 
dependence  in  the  far  infrared  is  good  in  the  3D  limit,  but  uncertain  in  2D.  The  control  required  of  p 
(~0.1  — 1%)  is  achievable. 
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Niklasson  and  Granqvist15  analyzed  old  data  of  Harris  et  a/.16-18  on  evaporated  Au  blacks  and 
concluded  that  anomalous  absorption  as  predicted  by  scaling  theory  is  observed  for  certain  samples  in 

the  mid-  to  far-infrared.  The  Au  blacks  are  tenuous  (p~0.002).  Since  the  particles  are  self-supporting 

0 

the  conducting  materials  must  be  close  to  pc.  The  particles  are  — 100  A  in  diameter.  The  thicknesses 
of  the  two  samples  of  interest  are  19  and  27  pm. 

NG  used  the  approximation 


^  -  1 


(12) 


to  obtain 


«  =  ^Re{cr(w;p)},  (13) 

where  9(w;p)  is  the  effective  conductivity.  Scaling  theory  predicts12  in  the  anomalous 

region.  Using  the  values  s=0.75  and  t  =  1.94  for  3D,  NG  calculated  the  exponent  0.72,  which  agrees 
well  with  the  exponents  0.703  and  0.717  obtained  by  converting  the  transmission  and  reflection  data  of 
Harris  et  al.  to  an  absorption  coefficient. 

Harris  and  Beasley27  state  that  the  approximation  given  by  Equation  12  does  not  apply  if 
A >  15  pm.  For  50  A  radius  Drude  Au  particles  in  the  anomalous  region,  |47ra/w|~|fl~19  and  188  at 
w  =  100  cm-1  for  3D  and  2D,  respectively.  However,  based  on  NG’s  Figure  1,  |4tt5/ui|~0.3  at  w  =  100 
cm-1.  Thus,  the  data  of  Harris,  et  al.,  as  analyzed  by  NG,  appear  to  be  consistent  with  the  use  of 
Equation  12.  We  conclude  that  there  is  an  inconsistency  between  the  prediction  based  on  the  Drude 
model,  which  accounts  for  the  FIR  properties  of  bulk  Au,  and  the  data  in  NG’s  Figure  1. 

The  correct  expression  for  the  frequency  dependence  of  a  for  Drude  Au  particles  is  given  by  Eq. 
8.  Using  the  same  values  of  s  and  t  used  by  NG,  we  obtain  exponents  of  0.86  and  3/4  for  3D  and  2D, 
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respectively.  Neither  exponent  provides  as  excellent  a  fit  to  the  data  as  does  0.72,  but  the  2D  exponent 
is  closer.  A  transition  from  3D  to  2D  percolation  behavior28-29  is  expected  as  pc  is  approached.  Using 
£~aAp-t,~L,  where  L  is  the  thickness  of  the  sample,  the  transition  is  estimated  to  take  place  at 
Ap~lxl0-4  for  L  =  19  /im.  It  is  unlikely  that  the  gold  black  deposits  are  so  close  to  pc.  Also,  the 
spectral  region  covered  by  the  data  does  not  fall  within  the  bounds  w*<w<w3  for  2D.  We  conclude 
that  the  scaling  theory  applied  to  Drude  Au  particles  is  not  able  to  explain  the  anomalous  frequency 
dependence  obtained  by  NG’s  analysis  of  the  data  of  Harris  et  al. 

Equation  8  provides  an  estimate  for  the  magnitude  of  a  very  close  to  pc.  For  3D,  a(w  =  100 
cm-1)  =  1.2xl0^  cm-1,  which  exceeds  the  measured  value  obtained  from  NG’s  Figure  1  by  a  factor  of 
seven. 

It  is  well-known  that  eddy  current  losses  (magnetic  dipole  absorption)  contribute  to  the  far 

1  30—31  8  32 

infrared  absorption  by  isolated  metal  particles  '  and  clusters  ’  .  The  question  of  whether  eddy 

current  losses  diverge  at  pc  has  been  addressed  for  both  normal  metal— insulator  and  superconducting 
composite  materials.  Both  numerical  work33-34  and  a  scaling  theory35  predict  divergent  response  at 
pc.  However,  a  recent  treatment36  that  takes  mutual  inductance  into  account  does  not  show  divergent 
behavior. 

The  scaling  theory  applied  in  this  paper  lead  to  a  homogeneous  effective  dielectric  function.  A 
scaling  theory  for  the  limit  k^^>l  based  on  geometrical  arguments3'  has  been  successfully  applied  to 
cermet  films  in  the  mid-  and  far-infrared38-42.  These  ideas  have  recently  been  developed  more  fully 
by  Yagil,  ef  a/.43  Rather  than  using  the  wavelength,  they  introduce  a  new  length  scale  L(u>),  the 
distance  a  carrier  diffuses  on  a  cluster  during  the  period  of  the  electromagnetic  wave  given  by  the 
theory  of  anomalous  diffusion44.  L(w)  decreases  with  increasing  frequency.  For  this  theory,  an 
effective  medium  treatment  is  not  valid  if  L(w)<C£,  and  the  optical  properties  of  the  sample  must  be 
calculated  by  summing  over  cells  of  size  L(w).  According  to  Yagil,  et  al.43,  L(w)  =  B£<,(A/£0)1/t2+^, 
where  ( =  £0|p  —  pc|' 1 ,  A  is  the  wavelength,  B  is  a  dimensionless  constant  of  order  unity,  and  0  =  1.5  in 
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three  dimensions  .  If  we  set  L(w)  =  £  with  £o=a=50  A  and  B  =  l,  we  obtain  the  crossover  frequency 
£?,.:=  i|p  —  pi:)ly^2+^^  =  (2.0xl04  cm_1)|Ap|3,1.  According  to  this  estimate,  a  description  based  on  finite- 

O 

size  scaling  is  required  to  calculate  the  far  infrared  properties  of  50  A  radius  Au  particles  near  pc 

We  conclude  that  dynamic  scaling  theory  predicts  interesting  frequency  and  composition 
dependence  for  the  far  infrared  absorption  coefficient  of  granular  metal-insulator  composite  materials. 
Experiments  on  cermet  films  and  other  composite  materials  in  the  far  infrared  with  composition  near 
pc  would  further  test  the  applicability  of  scaling  theories. 
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FIGURE  CAPTION 


Regions  of  behavior  of  the  far  infrared  absorption  coefficient  of  unsupported  50  A  radius  Au 
particles  near  pc  predicted  by  dynamic  scaling  theory.  The  horizontal  line  indicates  wr  =  l.  The  other 
solid  line  separates  the  region  of  anomalous  frequency  dependence  very  near  pc  from  the 

region  of  "normal"  frequency  dependence  (w<Cw*).  The  three  dashed  lines  are  estimates  of  the  limit  of 
the  region  of  applicability  of  the  quasistatic  approximation,  as  discussed  in  the  text,  a)  3D  b)  2D. 
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